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ABSTRACT 

Inertial  range  spectra  of  turbulent  flows  are  deduced 
from  the  equations  of  motion.   An  explanation  for  the 
validity  of  universal  equilibrium  theories  is  provided.   It 
is  shown  that  Fourier-space  cascades  and  physical-space 
intermlttence  are  complementary  aspects  of  the  same  reality. 
Spectra  of  the  form  E(k)  =  Ak~^  are  found,  where  k  is  the 
wave  number,  A  is  independent  of  k,  and  a  is  a  positive 
constant.   a  is  evaluated-,  for  the  Burgers  equation  a  =  2; 
for  the  Navier-Stokes  equations  in  both  two  and  three  space 
dimensions  a  =  2.5  ±  0.15- 
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Introduction 

The  velocity  field  of  a  fluid  in  turbulent  flow  is  the  sum 
of  many  components  of  widely  differing  scales.  A  fundamental 
objective  of  the  theory  of  turbulence  is  to  analyze  the  effect  of 
the  smaller  scales  on  the  larger  ones,  so  that  this  effect  can  be 
represented  by  an  appropriate  eddy  coefficient  in  a  closed  system 
of  equations  capable  of  solution  by  standard  means.   The  first 
step  towards  this  goal  is  the  study  of  the  smaller  scales  of 
motion  by  themselves;  in  many  cases  of  physical  interest  the 
important  scales  of  motion  belong  to  the  inertial  range,  which  we 
shall  now  define. 

Let  V  be  the  velocity  vector  in  a  turbulent  flow,  with 
components  v.  =  v.(x);  x, ,Xp,x^,  components  of  the  vector  x,    are 
the  space  variables.   Vector  quantities  will  be  underlined,  and 
averages  (whose  exact  meaning  remains  to  be  defined),  will  be 
denoted  by  bars  above  the  averaged  quantities. 

We  write 

v^(x)  =  v^(x)  +Uj^(x-)  . 
The  correlation  tensor  R.  .(Xj£)  is  defined  by  (see  [l]) 


f^li(2E.'Z)  =  u^(x)u  (x+r)  ; 


R.  ,  has  a  Fourier  transform  ^.  .  given  by 


ik.r 

.. ,  ,.  ..kU 

ij  -    J    ^iJ.- 


Rii(£)  -  /liifis)-'    ^^ 


where  the  dependence  of  R. .  and  ^. .  on  x  is  not  explicitly 
mentioned.   Let  k  be  the  length  of  the  vector  k,  and  let  dA(k)  be 
the  surface  element  on  a  sphere  S(k)  of  radius  k  centered  at  the 
origin.   Write 

$ij  =/    I^j(k)dA{k)  . 
S{k) 

The  energy  spectrum  E{k)  Is  defined  by 

E{k)  =  75-  tr.  •  (k)     (the  summation  convention 

IS  m  use) 

and  one  has  the  obvious  relations 

1  1 

Average  energy  at  x  =  2  u^(x)u^(x)  =  ^  R^^(O) 

00 
=  /  E(k)dk  . 
0 

It  can  be  shown  that  E{k)  Is  non-negative. 

We  assume  the  flow  Is  described  by  the  Navler-Stokes  equations 

^tZ.  ■•"  (Z.'Z^Z +SJ^a-d  p  =  vAv 
dlv  V  =  0 

where  t  Is  the  time,  p  Is  the  pressure  and  v  Is  the  viscosity. 

Prom  these  equations  It  can  be  seen  that  e,  the  rate  of  energy 

dissipation,  is 

00 


e  =  2v  /  k^E(k)dk 


2 
2k  E(k)  is  the  spectral  dissipation  function. 

In  flows  at  high  Reynolds  number  E(k)  differs  from  zero 

over  a  wide  range  of  k's,  and  furthermore  the  energy  spectrum  E(k) 

and  the  spectral  dissipation  function  2k  E(k)  are  approximately 

disjoint,  i.e.  there  exist  numbers  k.,  ,  k^  (not  uniquely  defined). 


such  that 


1  oo 


and 


/    E(k)dk  ^  /   E(k)dk  , 

0  'o 

oo  00 

1  k^E(k)dk  ^  I      k^E(k)dk 
k^  0 


^1  ^^  ^2 


The  intervals  (0,k-^),  (k-^jk^),  and  (k^jOo)  are  respectively  the 
energy  (or  energy-containing)  range,  the  inertial  range,  and  the 
dissipation  range. 

Attempting  to  find  general  expressions  for  E(x,k),  k  in  the 
energy  range,  is  both  Impossible  and  useless,  for  the  following 
reasons:   the  dependence  of  E  on  x  leads  to  problems  with  a  pro- 
hibitive number  of  independent  variables,  and  the  dependence  of  E 
on  the  initial  and  boundary  data  makes  the  existence  of  inter- 
esting general  results  unlikely.   The  study  of  E(k),  k  in  the 
energy  range,  should  be  carried  out  as  outlined  in  the  opening 
paragraph:   one  should  replace  the  effects  of  the  inertial  and 
dissipation  ranges  by  appropriate  eddy  coefficients  and  solve  the 
resulting  energy  range  equations  in. each  special  case.   Steps  in 
this  direction  have  been  described,  i.a.  in  [10]  and  [13]- 


The  study  of  the  inertial  and  dissipation  ranges  is  a 
sounder  undertaking.   Roughly  speaking,  energy  reaches  these 
ranges  after  a  number  of  non-linear  interactions;  the  influence 
of  the  boundary  and  initial  conditions  is  no  longer  appreciable; 
over  most  of  the  flow  the  spectrum  becomes  independent  of  x  and 
can  be  expected  to  embody  an  intrinsic  property  of  turbulence. 
It  is  of  course  a  task  of  a  theory  of  turbulence  to  provide  a  more 
credible  explanation  of  this  lack  of  dependence  between  the 
several  parts  of  the  spectrum.   When  the  Reynolds  number  is 
sufficiently  large  the  inertial  range  is  of  considerable  extent, 
and  its  influence  on  possible  eddy  coefficients  is  dominant. 
(See  e.g.  [11]  and  a  forthcoming  paper  by  the  present  author.) 
The  present  paper  is  devoted  to  the  study  of  the  inertial  range. 

The  Kolmogoroff  Law 

The  following  widely  quoted  result  in  the  theory  of  the 
inertial  range  is  due  to  Kolmogoroff.   Let  e  be  the  rate  of  energy 
dissipation, 

CO 


£  =  2v  /   k^E(k)dk  . 


0 

If  it  is  assumed  that  E(k),  k  In  the  inertial  range,  may  depend 
only  on  e  and  k  (and  not,  for  example,  on  the  viscosity  or  the 
Initial  data)  then,  on  dimensional  grounds. 


E(k)  =  constant,  e^/^  k"^'^^  . 


Theories  of  turbulence  tend  to  be  judged  by  their  ability  to 
reproduce  this  law  in  three-dir/  ■   .nal  flow;  yet  its  validity 
would  be  more  a  puzzle  than  an  explanation.   The  dimensional 
argument  is  independent  of  the  equations  of  motion  and  should  be 
valid  for  other  equations;  v;e  shall  see  in  the  next  section  that 
this  is  not  the  case.   The  validity  of  Kolmogorof f ' s  law  implies 
that  there  is  no  statistical  dependence  between  large  and  small 
scales  of  motion  and  that  at  high  Reynolds  numbers  the  rate  of 
energy  dissipation  is  Independent  of  viscosity  (see  the  excellent 
discussion  in  [12]).   It  is  far  from  obvious  that  this  should  be 
true.   The  hypotheses  on  v/hich  Kolmogorof  f '  s  lav;  is  based  are  in 
no   way  necessary  for  the  existence  of  an  Inertial  range.   It  is 
therefore  of  interest  to  evaluate  the  inertial  range  spectr-um  from 
a  dynamic  theory,  i.e.  using  the  I^avier-Stokes  equations  or  a 
credible  theory  of  turbulence. 

Kraichnan  [7]  has  performed  such  a  computation  on  the  basis 
of  his  direct  interaction  approximation  and  obtained  for  three- 
dimensional  flov; 

E(k)  =  constant -k"^^ 

Kraichnan' s  later  work  is  not  relevant  in  this  context,  slnct  it 
is  based  on  the  assumption  that  Kolmogorof f ' s  law  is  valid  and 
does  not  afford  an  independent  check  of  the  latter.   Tlie  author 
knows  of  no   other  previous  attempt  to  verify  Kolmogorof  f  s  lav/. 
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Several  inertlal  range  spectra  for  two-dimensional  turbu- 
lence have  been  proposed  of  late,  also  on  the  basis  of  dimensional 
arguments.   The  energy  spectrum 


E(k)  =  constant- k 


has  been  derived  by  Leith  and  by  Kraichnan  (see  [8],  [9]).   It  is 
interesting  to  note  that  all  proposed  laws  are  of  the  form 


E(k)  =  constant'k    ,    P  constant 


A  discussion  of  these  expressions  will  be  presented  later. 

The  method  described   in  this  paper  can  be  summarized  as 
follows:   The  Fourier  transform  u(k,t)  of  u(x,t)  is  taken.   (An 
outline  of  a  justification  of  this  step  can  be  found  in  [  1  ].) 
An  equation  of  motion  for  u(k,t)  is  derived  from  the  Navier-Stokes 
equations  and  the  assumption  that  |k|  is  the  inertlal  range.   The 
resulting  equations  are  essentially  transforms  of  the  Euler  equa- 
tions, subject  to  certain  restrictions  on  the  possible  form  of  the 
solution.   Stationary  solutions  of  these  equations  are  then  sought. 
The  reason  for  singling  out  stationary  solutions  is  best  brought 
out  by  an  analogy  with  a  finite  mechanical  system:   Consider  such 
a  mechanical  system  with  N  degrees  of  freedom  denoted  by  the 
variables  a,  ,  .  .  .  ,cx  ,  and  whose  motion  is  described  by  equations  of 
the  form 

da. 

(1)  ^  =  F.(a^,...,a^) 


with  initial  conditions 

a^(0)  =  g^ 

where  the  g^  are  random  with  an  unknown  distribution  law.   If 
a  stable  stationary  solution  a.  of  (i)  can  be  found,  to  which  all 
solutions  will  tend,  then  at  time  t  large  enough  the  average 
values  of  a.  will  be  approximated  by  a. ,   The  appropriate  average 
is  obviously  an  average  over  an  ensemble  of  realizations  of  the 
system. 

For  the  sake  of  clarity,  we  shall  first  apply  our  method  to 
the  one-dimensional  Burgers  equation.   This  is  the  object  of  the 
following  section. 

The  Burgers  Equation  and  its  Inertlal  Range 

Consider  the  Burgers  equation 

1      2    12 

(2)       S.u+TT-  5(u) Su  =  0,    R  constant  . 

t    2   X       R   X       ' 

The  behavior  of  the  solution  of  (2)  can  be. studied  by  fairly 
standard  methods.   (In  fact  (2)  can  be  reduced  to  a  linear  equation 
and  solvea,  but  this  will  be  of  little  use  in  the  present  context.) 
Solutions  of  the  "invlscld"  equation 

S^u  +  I  ^^[^^)    =   0 

are  obtained,  and  rendered  single-valued  by  apt  Insertion  of  shocks. 
The  structure  of  these  shocks  can  be  determined  and  the  energy 
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spectrum  as  well  as  the  dissipation  function  approximately 
evaluated  (see  [12]).   An  inertlal  range  is  found,  with 

(3)  E(k)  =  Ak"^ 

where  A  is  independent  of  k  (but  does  depend  on  the  structure  of 
the  initial  data).   The  Kolmogoroff  hypothesis  is  thus  invalid. 

We  now  proceed  to  rederive  (3)  by  a  method  capable  of  exten- 
sion to  fluid  flow  in  several  space  dimensions.  Let  u(k,t)  be  the 
spatial  Fourier  transform  of  u(x,t).   It  satisfies 

(4)  S^a(k,t)  =  -  ^y  u(k',t)u(k-k',t)dk'  -  J  k^u(k,t)  . 

The  real  character  of  u  is  reflected  in  the  condition 

(5)  0(k,t)  =  u*(-k,t) 

where  the  asterisk  denotes  a  complex  conjugate  function.   Let  k 
belong  to  the  inertial  range.   By  hypothesis,  the  dynamic  effect 
of  the  viscosity  R~   on  u(k,t)  is  negligible  compared  to  the  inertia 
terms.   The  last  term  in  { 4  )  can  be  neglected,  yielding 

(6)  B^u  =  -  -^y^u(k',t)u(k-k',t)dk'  . 

It  should  be  emphasized  that  the  neglect  of  the  dissipation  term  in 
(  ^  )  is  essentially  different  from  the  Kolmogoroff  assumption  that 

E(k)  is  independent  of  R~  .   The  equations  are  strongly  non- 
linear, and  a  strong  dependence  on  R  in  other  parts  of  the  spectrum 
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may  strongly  affect  the  amplitude  of  E  in  the  inertlal  range.   We 
shall  in  fact  find  spectra  different  from  Kolmogorof f ' s.   The 
fact  that  energy  is  dissipated  at  large  k  leads  one  to  expect  that 


lim  |u(k) I  =  0 
k— »•  00 


This  expectation  will  be  fulfilled. 

We  note  that  equation  (  6  )  is  invariant  under  the  transforma- 
tion 


^  /S 


{  7) 


u'  =  u/U 


k'  =  k/K 


t'  =  t(UK)  . 

In  particular,  choose  K  to  be  a  wave  number  in  the  inertial  range, 
and  U  =  u(K).   Having  assumed  that  the  energy  is  contained  mostly 
in  the  energy  range,  we  expect  that  u/U  will  become  large  as  k/K 
decreases;  if  (  6  )  is  considered  to  be  an  equation  for  u',  this 
leads  to  the  condition 


(  8) 


lim  I  u(k)  I  =  oo 
k-^0 


The  integral  in  (  6  )  is  necessarily  singular,  and  must  be  inter- 
preted in  the  obvious  generalized  sense 


u(k' )u(k-k' )dk'  =  lim 


-e 


k-e 


+ 


+ 


oo 


k-fe 


u(k')u(k-k' )dk' 


One  can  write  (  6  )  and  (  8  )  in  the  symbolic  form 


(  9  )  5^a  =  F[Q] 


where  F[u]  =  F[u](k)  is  a  functional  of  u  and  a  function  of  k, 
We  now  seek  steady  solutions  of  (  9  ),  i.e.  solutions  of 


■ik  /  u(k',t)u(k-k' )dk'  =  0 


In  the  present  one-dimensional  case  such  a  solution  is  readily 
found.   Let  u(x,t)   be  the  function  of  which  this  solution  is  a 
transform.   Obviously,  we  have 


5  (u^)  =  0 
x^   ^ 


and 


«v 


(10)  u(x)  =  ±C  ,    C  constant  , 

different  signs  being  assumed  on  separate  intervals  of  the  x-axls. 
From  this  one  concludes  that 

( 11)  u(k)  =  constant'k" 

and 

E(k)  -  Ak'^  . 

The  solutions  ( 10)  and  ( 11)  are  worth  examining,   u  represents 
an  array  of  shocks  with  arbitrary  locations,  i.e.  an  intermittent 
flow.   u  describes  a  wave  number  space  cascade;  Indeed,  if 
0  <  k'  <  k,then 

10 


1    1   >  0 


k'  k-k' 
and  if  k'  >  k  or  k'  <  0 

1    1 


k'  k-k 


r  ^  0 


Thus  u(k,t)  is  fed  by  modes  with  wave  numbers  smaller  than  k  and 
depleted  by  interactions  with  wave  numbers  larger  than  k. 
Intermittent  flows  and  wave -number  space  cascades  appear  as 
complementary  aspects  of  the  same  reality,  not  unlike  waves  and 
particles  in  quantum  theory. 

A  Numerical  Method 

It  will  now  be  shown  how  (11)  can  be  derived  by  a  method 
applicable  to  multidimensional  problems  as  well.   We  seek 
stationary  solutions  of  (9);  these  can  be  determined  only  up 
to  an  arbitrary  constant.   If  the  energy  spectrum  E{k)  is  to  be 
unique  except  for  a  multiplicative  constant,  we  must  have 

E{k/K)  =  u(k/K)u*(k/K)  =  constant  •u(k)u''{k)  =  constant  •E(k)  . 

Put 

u(k)  =  a(k)  exp  (i(|)(k)),    a{k)  real  and  positive  . 

It  follows  that 

a(k/K)  =  constant .a(k) 

i.e.,  in  view  of  (8 ) 

(12)  a(k)  =  constant. I k| "^  ,    P  >  0  . 
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We  are  in  effect  looking  for  similarity  solutions  of  (6  ). 

Consider  now  the  possible  forms  of  ij)(k).  ^   must  be  an  odd 
function  of  k,  to  satisfy  (  5  ).   For  F[u]  to  be  finite, 
(t)(k' )  +  (j)(k-k '  )  must  change  by  tt  at  k'  =0  and  k'  =  k;  this  will 
cause  the  singularities  in  u(k')u{k-k')  to  cancel.   Furthermore, 
(t){k)  must  be  such  that 

(13)  ^(k')  +(j)(k-k')  =  g{k)  +h(k'/k) 

where  g  depends  only  on  k  and  h  only  on  k'/k.  The  significance 
of  this  last  condition  will  become  apparent  shortly.  Functions 
(()  which  satisfy  all  these  conditions  are 

(14)  (t)(k)  =  ak  -sgn(k)ir/2  ,     a  an  arbitrary  constant, 

sgn{k)  =  k/|k|  . 

This  represents  a  natural  choice  of  (j)  in  view  of  the  following 
observation:   Consider  an  initial  value  problem  for  equation  (4  ) 
with  Initial  values  consisting  of  a  single  Fourier  mode, 

u(k)  =  5(k-l)  exp  {li/   sgn(k)) 
i.e. 

u  =  cos  {x+Tp)    . 
It  can  be  verified  that  for  all  t 
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u(k,t)  =  ')         b(k,t)6(k-n)  exp  {lrn^+(n-l)  sgn{n)T7-/2) 


n 


)    b(n,t)6(k-n)  exp  ( ia^  -  sgn(n)7r/2)  , 
n 

a  =  ^  +  7r/2  . 


This  form  of  ^   has  already  been  discussed  in  [4  ]. 

Given  ^,    there  remains  only  the  task  of  determining  p.   If 
F[u](k  )  =  0  for  some  k  ,  then  F[u](k)  =  0  for  all  k;  this  is 
ensured  by  the  condition  (  13)  and  follows  from  the  following 
identities  : 


F[u](k)  =  constant-/  Ik'T^jk-k 


-P 


.  exp  (iak  -  i(7r/2)  (sgn(k'  )  +  sgn(k-k' ) )  )dk'- 

-e   k-e   00  -1 


=  exp  (iak)  lim 

e-*.0 


+ 

•00    £ 


k+e 


(Ik'r^lk-k'r^)dk'  . 


Putting 


~  k  ^  ' 


we  obtain 


1    k   -28 
F[u](k)  =  exp  (iak)  ^  (  °)   ^  lim 

o        e-»-0 


•  ek  /k   k  -ek  /k   oo 
o'     p  o    o'    n 


+ 


-00 


ek  /k 

o^ 


k  +ek'  /k 
o   o' 


or 


k'T^lk^-k'T^dk" 


k   -26-1 
F[u](k)  =   exp  (ia(k-k^))(^)     F[a](k^) 


15 


The  problem  is  therefore  reduced  to  finding  p  such  that 


(15)  F[u^](l)  =  0 


where 

u^  =  k/|k|^+^  =  e^^/^|k|-P  exp  {-sgn(k)Tr/2)  ,  ' 

and  a  was  set  equal  to  zero  without  loss  of  generality.   The 
point  k  =  1  is  chosen  so  that  all  the  u  (k  )  have  the  same 
amplitude,  and  information  about  the  location  of  the  zero  can  be 
obtained  from  the  variation  of  the  residuals 


rp  =  la^G^d)  =  iP[a^](i) 


In  Table  I  the  results  of  such  a  computation  are  exhibited.   The 
integration  was  performed  as  follows:   Let 

f(k')  =  uP(k')u^{l-k')  ; 
the  generalized  Integral 

+00 

I        f(k')dk' 

-OD 

is  replaced  by  the  equivalent  expression 
1 


+00 

/   g(k')dk'  ,    g{k')  =  Y^     f{j+k') 

0  J=-°° 

g  can  be  shown  to  be  a  bounded  function  of  k'.   g{k')  is  approxi- 
mated by 


Ik 


Table   I 


Minimization   of  F[u'^] 


M   =    100;    e    =    10    -^°;    k      =    1 
'  '      o 


p 

'^P,N 

''P,N 

0.90 

.9747 

.096 

0.95 

.479^ 

.023 

0.98 

.1896 

.0036 

0.99 

.0946 

.000092 

1.00 

-.0044 

.0000020 

1.01 

-.1027 

.0010 

1.02 

-.1984 

.0040 

1.05 

-.5056 

.026 

1.10 

-.9955 

.10 

15 


+N 
^^        -N 

where  N  is  chosen  so  that 

|f (N-1) I  +  |f (-N+1) I  ^  e  ,    e  small  . 

g„  was  then  integrated  by  a  Monte-Carlo  procedure.   In  Table  I 
rp,    is  the  result  of  this  integration  and  Co  ^  the  corresponding 
variance.   The  results  show  that  p  =  1  is  the  appropriate  value 
of  P  and  the  result  of  the  preceding  section  is  recovered.   We 
have 

E(k)  =  Ak"^  . 

The  question  of  the  stability  of  these  stationary  solutions 
remains  to  be  answered.   The  conclusion  that  they  are  stable  can 
be  drawn  from  the  existence  of  known  solutions  of  (  4  )  which  for 
large  t  and  large  k  have  the  spectrum  (  3  )j  see  [  2  ],  [  3]  and  [12]. 
The  numerical  evidence  drawn  from  finite  difference  integrations 
of  (  4)  also  points  to  this  conclusion.   The  phase  function  ^   appears 
clearly;  in  fact  this  is  how  it  was  originally  found.   With  good 
will  and  a  large  number  of  Fourier  modes  in  the  calculation,  the 
spectrum  (  3  )  can  also  be  fleetingly  observed;  it  is  soon  destroyed 
by  the  combined  effects  of  a  finite  number  of  Fourier  modes  and  a 
finite  time  step. 

The  reason  for  the  stability  of  these  solutions  is  easy  to 
visualize  in  x-space.   An  unsteady  solution  will  consist  of  a 
number  of  shock  waves  of  differing  intensity  and  differing  speeds. 
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Some  will  coalesce,  some  will  disappear  to  infinity,  until  there 
remain  only  shocks  of  equal  intensity  and  speed. 

The  Inertial  Range  in  Turbulent  Incompressible  Flow 

Let  u(k, t),  with  components  u  (k,t),  be  the  Fourier  trans- 
form of  the  fluctuating  velocity  vector  u(x,t);  the  Fourier  trans- 
form of  the  Navier-Stokes  equations  is,  in  non-dimensional  form, 

(16)    \n^(]i.t)    =  -lk,P„^ag^  -  ^  "a  '    "a'-"^'  =  <"^' 


where 


Q   =  /  u  {k-k')G  (k' )dk'  , 

k  k 

p   =5 ^^     (5    the  Kronecker  delta)  , 

a/    a/    ,2       ^    ay 


and  the  summation  convention  is  in  force.   The  pressure  has  been 
eliminated  through  the  use  of  the  equation  of  continuity 


(17)  k  u  =  0  , 


giving  rise,  in  the  well-known  manner,  to  the  projection  P^  . 
R  is  the  Reynolds  number.   We  shall  now  study  inertial  range  solu- 
tions of  these  equations.   As  before,  the  viscous  dissipation  term 
(-k^/R)u  is  dropped,  and  the  invariance  of  the  equation  under  the 
transformation 


17 


u'  ==  u/U 
(18)  k'  =  k/K 

t'  =  tUK 

is  noticed;  this  leads  to  the  conditions 


(19)      lim  |Q  (k)|  =00  ,     l^^(lS)l  bounded  for  |k|  ^  0    . 
|k|-*0  "  " 


As  before,  the  integrals  Q  o  must  be  interpreted  in  a  generalized 
sense.   We  rewrite  (16)  in  the  symbolic  form 

(20)  B^u  =  F[u]  =  F[0](k) 

and  seek  steady  solutions  of  this  equation.   Consider  the  case  of 
two-dimensional  flow.   The  equation  of  continuity  (l?)  is  satisfied 
if 

^1  ""  -k2p(k) 

Ug  =  k-j^p(k)  . 

p  is  a  stream  function.   From  (19)  we  obtain 

(21)  lim  |p(k)  I  =  00  ,    k  =  |k|  . 

k— 0 

A  function  p(k)  invariant  under  the  transformation  (18),  satisfying 
(19),  and  ensuring  the  boundedness  of  F[u],  is 


P^iS)  =  HlkT^ 


iLk 
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where  Lk  is  an  arbitrary  linear  function  of  k  and  H(k)  is  a 
homogeneous  polynomial  in  k..  and  kp  which  vanishes  only  when  k  =  0. 
If  H(k)  is  such  that  F[u]  vanishes  at  all  points  of  a  circle  in 
the  (k-,  jko)  plane,  then  F[u]  vanishes  everywhere  by  construction. 
F[u]  can  vanish  only  if  H(k)  reflects  the  fact  that  equation  (l6) 
Is  invariant  under  the  interchange  of  the  k-,   and  kp  axes; 
nevertheless  H(k)  contains  an  infinite  number  of  parameters,  to  be 
determined  by  satisfying  the  infinite  number  of  equations 

(22)  F[u](k)  =  0  ,    |k|  =  constant  . 

This  will  be  done  as  follows:   H(k)  will  be  chosen  to  be  invariant 
under  interchange  of  k-,  and  kp  and  under  the  transformation 

K  =  -^1 

kg  =  -kg  . 

It  is  then  sufficient  to  satisfy  (  22)  on  the  portion  of  the  unit 
circle  lying  between  the  lines  k^  =  0  and  k^-kp  =  0. 
On  the  unit  circle  all  the  functions  H(k)  have  comparable 
magnitude.   Our  effort  will  thus  be  directed  towards  the  minimiza- 
tion of  appropriate  residuals  of  equation  (22). 

Two  Dimensional  Flow 

We  now  proceed  to  carry  out  the  minimization  of  P  in  the  case 
of  two  dimensional  flow.   We  choose  n  equidistant  points 
_s^  =  (cos  e.   ,sin  9.) ,    i  =  1,.  .  .  ,n  on  the  unit  circle,  6^=  6^/2  +  19^, 
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0  <   6 .    ;^  ir/k;    If  F  can  be  made  to  vanish  on  those  n  points.  It 
will  vanish  by  symmetry  on  8n+4  points  of  the  whole  unit  circle. 
We  define  the  total  residual  relative  to  u  to  be 

Total  residual  =^11   UPit^.^  (  s.^  )  I  +  1^2^-^  ^-^  ^  '  •*  ' 

i 

The  total  residual  is  minimized  first  within  a  one  parameter  family 
of  functions  H,  then  within  families  of  functions  with  a  larger 
number  of  parameters.   Each  new  parameter  is  introduced  in  a  way 
which  makes  use  of  the  information  previously  obtained. 

The  point  9  =  0  is  to  be  avoided  on  computational  grounds: 
at  this  point  the  singularities  of  Q  o  coalesce  and  an  added 
numerical  effort  would  be  required.   At  the  point  9   =  ir/k,   F  is 
always  zero  by  symmetry.   Write 


f^   g(k')  =  Ug(k' )u  (s-k' )  ,   _s  =  (cos  9, sin  ^9)  ; 


the  generalized  integral 

is  replaced  by  the  equivalent  expression 

where 

D  =  Ck-[,k2  I  0  <  k^  <  cos  9,    0  <  k^  ^  sin  9) 
and 
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+CX)     +00 


§57  e   =  ^ZI   ZZ  ^5^  0(i  cos  e  +  k'   j  sin  0+k')  . 
^  '     i  =  -oo  j  =  -oo    ^ '  ^ 

is  now  a  bounded  function  of  k'.   The  infinite  sum  is 


^67,0 
approximated  by 

^bj,en  =  .IZ^  ^67,0 

where  N  is  determined  by 

>     (  I  f-i-i  ,a|  +  I  f-|o  fli  +  I  foQ^-fll  )  1  e  ,    e  small. 
l+j=N  ^  '       "^ 

It  is  convenient  to  choose  e  to  be  a  multiple  of  V  =  cos  0   sin  0, 
the  area  of  D. 

The  integration  of  g  has  been  carried  out  by  a  straight- 
forward second  order  integration  formula,  with  M  points  in  each 
direction. 

The  first,  one  parameter,  family  of  functions  H  was  chosen 
to  be 

H(k)  =  k  ^  ,    k  =  /  ^i  +^2   "  '-'  • 

The  reason  for  this  choice  is  not  only  the  analogy  with  the  Burgers 

equation,  but  also  the  following  argument:   for  d,u  to  vanish,  the 

vector  with  components  kc-Qj,   must  be  in  the  null-space  of  the  operator 

P   ,  i.e.  it  must  be  orthogonal  to  the  unit  circle.   The  closer  u 
ay'  °  — 

comes  to  rotational  invariance,  the  more  likely  it  is  that  this 
orthogonality  will  hold.   There  exist  no  non-zero  rotation-invariant 
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vectors,  but  in  an  intuitive  sense  this  choice  of  H  brings  one 
close  to  rotational  invariance.   Additional  parameters  are  intro- 
duced as  indicated.   The  spectrum  is  given  by 

77" 

E(k)  =  constant-k  /   u(k)u*(k)d0  ,    cos  0  =  k^/|k|  ,    |k|  =  k  . 

0 

One  may  notice  that,  as  in  the  one-dimensional  case,  the 
solutions  we  shall  obtain  represent  a  flow  of  energy  from  low  to 
high  wave  numbers.   They  are  the  Fourier  transforms  of  distribu- 
tions of  vortices. 

The  functions  gj-   g  are  not  well  behaved,  in  particular  they 
are  subject  to  large  oscillations  around  the  points  ]£'  =  0  and 
k'  =  _s.   Numerical  integration  would  therefore  appear  to  be  a 
difficult  undertaking.   Fortunately  the  functions  g(k  '  )  +  g(  s_-k' ) 
are  reasonably  smooth,  and  reliable  results  can  therefore  be 
obtained  by  an  integration  formula  which  is  symmetric,  i.e.  invari- 
ant under  the  transformation  g(k')  -*   g(s_-]£' )  for  each  _s.   e  must  be 
so  small  that  the  effect  of  further  change  in  e  on  the  location  of 
the  smallest  residual  is  small;  more  stringent  requirements  cannot 
be  met  in  reasonable  computing  time.   Tables  II  and  III  demonstrate 
that  the  effects  of  a  change  in  the  number  of  points  M  in  the  inte- 
gration formula  and  of  a  change  in  e  are  not  major.   It  can  also  be 
verified  that  the  effect  of  a  change  in  the  number  of  points  n  on 
the  circle  on  which  the  residual  is  computed  is  minor  when  n  >  3' 

The  residuals  are  small  compared  to  what  they  are  in  the  case 
of  the  Burgers  equation.   One  cannot  conclude  from  this  that  all 
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Table  II 


Minimization  of  P_[u] 

2^1 
H  =  k    ;  e  =  O.O5V;  M  =  5;  n  =  3 


Total  residue 


^1 

xv^i^CLJ.   J.  c:  0 -L  VJ.  U.  C 

u 

1.425 

.0180 

2.7 

1.400 

.0099 

2.6 

1.375 

.0059 

2.5 

1.350 

.0048 

2.4 

1.325 

.0089 

2.3 

1.300 

.0147 

2.2 

1.275 

.0206 

2.1 

1.250 

.0257 

2.0 

1.225 

.0324 

1.9 

1.200 

.0328 

1.8 

1.175 

.0462 

1.7 

1.150 

.0557 

1.6 

1.125 

.0675 

1.5 
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Table  III 


Minimization  of  F[u] 
H  =  k    ;  e  =  0.02V;  M  =  6;  n  =  3 


P-      Total  residue  a 

1.550  .0678  5.2 

1.525  .05^^  5.1 

1.500  .0435  3.0 

1.475  .0338  2.9 

1.450  .0256  2.8 

1.425  .0187  2.7 

1.400  .0126  2.6 

1.375  .0068  2.5 

1.350  .00332  2.4 

1.325  .00330  2.3 

1.300  .0069  2.2 

1.275  .0112  2.1 

1.250  .0156  2.0 

1.225  .0205  1.9 

1.200  .0258  1.8 

1.175  .03191  1.7 


24 


our  trial  solutions  are  in  some  sense  nearly  solutions  of 
FLuJ  =  0;  it  should  be  remembered  that  if  we  change  the  circle  on 
which  the  minimization  is  carried  out  some  of  the  residuals  may 
be  considerably  amplified.   This  remark  may  be  used  for  computing 
upper  bounds  and  lower  bounds  for  the  exponent  a  of  the  inertial 
range  spectrum. 

In  the  numerical  tables  a  is  the  inertial  range  exponent 
associated  with  the  given  trial  function  H,  i.e.  the  inertial 
range  spectrum  would  be  Ak"*^  if  the  trial  function  H  under  con- 
sideration annihilated  F. 

Table  IV  demonstrates  that  the  addition  of  parameters  to  H, 
i.e.  the  use  of  more  complicated  trial  functions  H,  causes  only 
marginal  changes  in  a.   The  conclusion  from  the  tables  is  that 

a  =  2.3  ±  0.15  . 

It  is  hoped  that  future  work  will  yield  a  more  accurate  determina- 
tion.  The  value  of  a  is  almost  exactly  in  the  middle  between  the 
value  a  =  3  found  by  Leith  [9]  and  Kraichnan  [8]  and  the  value 
(J  =  5/3  found  by  Kolmogoroff.   Both  Kraichnan  and  Leith  state  that 
the  values  a  =  5/3  and  a  =  3  can  both  apply  to  two-dimensional 
turbulence,  and  therefore  to  a  limited  extent  the  present  calcula- 
tion confirms  their  results. 

The  coefficient  A  in  the  Inertial  range  spectrum  depends  on 
the  nature  of  the  large  scale  flow  and  cannot  be  determined  by  the 
present  theory.   In  particular  one  can  provide  examples  of  large 
scale  flows  which  do  not  give  rise  to  an  inertial  range.   One 
example  is  the  periodic  flow 
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Table  IV 
Minimization  of  F[u] 


2    2,^1  .   .2Pn   .   .2P. 


H(k)  =  (k^  +k2)  ^(|k^|   1+  Ikgl   2);  e  =  O.OIV;  M  -  5;  n  =  4 


Pl 

^2 

total 
residue 

0 

1.40 

0. 

.0126 

2.6 

1.39 

,01 

.0180 

2.6 

1.38 

.02 

.0476 

2.6 

1.3250  0.  .0076  2.5 

1.3150     .01  .0236  2.3 

1.3125  0.  .0052  2.25 

1.3025     .01  .0258  2.25 

1.3100  -.01  .0307  2.20 

1.3000  0.  .0039  2.20 

1.2900       .01  .0277  2.20 

1.2850  -.01  .0240  2.10 

1.2750  0.  .0044  2.10 

1.2650  .01  .0323  2.10 

1.2600  -.01  .0172  2.00 

1.2500  0.  .0092  2.00 

2.2400  .01  .0366  2.00 
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(23)     u  =  (C  cos  ax,  sin  aXp, -C  sin  ax^  cos  aXp )  , 

C,  a,    constant  . 

-t/R 
u  is  a  solution  of  Euler's  equations,  ue   ^   is  a  solution  of  the 

Navier-Stokes  equations  for  all  time  t.   Energy  is  not  transmitted 
to  the  higher  modes  and  therefore  A  =  0. 

Three-Dimensional  Flow 

It  can  be  seen  that  a  two-dimensional  solution  of  the  equa- 
tion F[u.]  =  0  is  also  a  solution  of  this  equation  in  three- 
dimensional  space;  this  solution  satisfies  the  condition 

11m  I  u(  k)  I  =  00 
k-fO 

and  for  k  7^  0  has  a  finite  spectrum  E(k);  these  conditions  are 
natural  generalizations  of  conditions  (19) >  fully  consistent  with 
the  assumption  that  k  lies  in  the  inertial  range. 

An  attempt  was  made  to  find  other  solutions  of  F[u]  =  0  in 
three  dimensions  by  an  appropriate  sequence  of  minimizations.   None 
were  found.   The  reason  for  the  absence  of  genuinely  three- 
dimensional  solutions  can  be  explained  as  follows:   because  of  the 
equation  of  continuity  (17),  u  must  be  of  the  form 

u  =  k  X  p 

where  x  denotes  a  vector  product  and  p  is  a  vector  potential.   For 
any  choice  of  p,  u  lies  in  a  plane  orthogonal  to  p,  and  so  does  the 
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vector  w  with  components  w  =  ■^5'^6-y'  ^^   ^^^   ^^  seen  that  w  7^  0 
for  any  reasonable  choice  of  £;  thus  for  F(^[u]  =  ^^b^ay'^by  ^° 
vanish  w  must  be  in  the  null  space  of  the  projection  P   ,  i.e. 
must  be  orthogonal  to  the  unit  sphere.   This  is  impossible. 

In  physical  terms  the  fact  that  the  solutions  are  two- 
dimensional  can  be  explained  as  follows:   The  inertial  range  flow 
consists  of  a  random  collection  of  vortices;  in  the  neighborhood 
of  each  vortex  the  fluid  flows  in  a  plane  orthogonal  to  the  vortex, 
i.e.  it  is  two-dimensional. 

There  are  clearly  differences  between  the  two-dimensional 
and  three-dimensional  cases  which  account  for  the  disorderly 
appearance  of  three-dimensional  flow.   The  orientations  of  the 
vortices  (the  normals  to  the  planes  of  the  flows)  are  random,  and 
it  can  be  argued  that  the  amplitudes  of  the  motions  are  larger  in 
three  dimensions.   In  particular,  no  three-dimensional  flow  similar 
to  (23)  is  known,  such  that  the  amplitude  of  the  corresponding 
inertial  range  flow  Is  zero;  I  conjecture  that  none  exists. 

The  fact  that  in  the  small  the  solution   of  the  equations  of 
motion  is  two-dimensional  should  be  less  surprising  since  we 
already  know  of  a  three-dimensional  flow  which  left  to  itself  be- 
comes two-dimensional;  this  happens  in  a  convection  problem,  see 

e.g.  [5]. 

The  inertial  range  spectrum  of  three-dimensional  turbulence 

is  thus  of  the  form 


E(k)  =  Ak  ^  ,    a  =  2.3  ±  0.15  . 
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This  result  is  at  variance  with  Kolmogorof f ' s  law,  which  Is  thus 
found  to  be  Invalid  In  all  cases. 

The  lack  of  validity  of  Kolmogorof f ' s  law  is  not  very 
surprising.   It  is  not  necessary  to  attach  undue  Importance  to 
results  obtained  by  dimensional  analysis,  since  the  possible  rela- 
tions between  the  various  properties  of  Inertlal  range  flow  are 
only  dimly  understood.   It  should  be  remembered  that  in  two  space 
dimensions  dimensional  arguments  yield  two  distinct  spectra. 

Inertlal  range  flow  should  be  visualized  as  follows:   vortices 
are  constantly  being  shed  from  regions  of  shear  in  the  main  flow. 
Rapidly  they  assume  an  asymptotic  form  with  the  spectrum  we  have 
just  computed.   Their  amplitude  decays  under  the  influence  of 
viscosity,  and  they  are  constantly  replaced  by  new  vortices.   This 
picture  provides  an  explanation  of  the  independence  of  the  form  of 
the  inertlal  range  spectrum  from  the  form  of  the  large-scale  motion. 

It  may  be  worth  pointing  out  that  equation  (20)  is  derived 
from  equation  (l6)  through  an  argument  strongly  reminiscent  of 
arguments  used  in  boundary  layer  theory;  in  fact  equation  (20)  can 
be  obtained  by  expanding  u  formally  in  Inverse  powers  of  a  Reynolds 
number  R    characteristic  of  inertlal  range  flow.   The  search  for 
similarity  solutions  is  natural  In  view  of  the  absence  of  well- 
defined  length  scales  in  inertlal-range  flow,  and  is  of  course  a 
common  device  in  boundary  layer  theory. 
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Comments 

The  present  paper  has  arisen  from  an  attempt  to  solve 
numerically  the  Green-Taylor  problem  [l4].  I.e.  from  an  attempt 
to  study  properties  of  turbulence  by  examining  numerical  solutions 
of  Initial  value  problems  for  the  Navler-Stokes  equations.   Such 
attempts  are  doomed  to  failure.   The  number  of  mesh  points 
required  is  huge  (an  estimate  can  be  obtained  from  the  fact  that 
the  sum  defining  g^  q^t   contains  thousands  of  terms),  and  the 
accumulated  errors  in  the  Integration  in  time  are  unavoidably  large 
(see  [6]).   In  fact,  the  inertial  range  spectrum  obtained  by  such 
means  is  devoid  of  physical  significance  and  depends  mainly  on  the 
difference  scheme  employed.   The  study  presented  in  this  paper 
avoids  these  difficulties.   It  does  however  verify  a  fundamental 
premise  of  numerical  work  on  the  Green-Taylor  problem:   it  shows 
that  the  structure  of  turbulence  can  be  studied  by  examining 
properties  of  a  few  typical  flows. 

The  premises  of  the  present  work  are  however  at  variance 
with  the  assumptions  of  most  statistical  theories  of  turbulence. 
The  only  statistical  argument  used  is  the  assertion  that  if  a 
number  of  solutions  tend  to  a  common  limit,  so  does  their  average. 

Finally,  the  present  paper  yields  not  only  an  approximation 
to  the  form  of  the  inertial  range  spectrum,  but  also  further 
information  about  the  nature  and  structure  of  inertial  range  flows; 
this  can  be  used  for  evaluating  eddy  coefficients.   The  manner  in 
which  this  is  to  be  done  will  be  the  subject  of  a  forthcoming 
paper. 
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Appendix 

For  the  convenience  of  anyone  who  would  wish  either  to 
check  or  to  extend  the  calculations  presented  In  this  paper, 
the  program  used  to  obtain  the  two  dimensional  results  Is 
reproduced.   No  particular  effort  has  been  expended  to  make 
the  program  very  elegant  or  very  efficient. 
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